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We report on a combined scanning tunneling microscopy (STM) and density lunctional theory 
(DFT) based investigation ol Co atoms on Ge(lll)2xl surfaces. When deposited on cold surfaces, 
individual Co atoms have a limited diffusivity on the atomically flat areas and apparently reside on 
top of the upper 7r-bonded chain rows exclusively. Voltage-dependent STM imaging reveals a highly 
anisotropic electronic perturbation of the Ge surface surrounding these Co atoms and pronounced 
one-dimensional confinement along the 7r-bonded chains. DFT calculations reveal that the individual 
Co atoms are in fact embedded in the Ge surface, where they occupy a quasi-stationary position 
within the big 7-member Ge ring in between the 3 rd and 4 th atomic Ge layer. The energy needed for 
the Co atoms to overcome the potential barrier for penetration in the Ge surface is provided by the 
kinetic energy resulting from the deposition process. DFT calculations further demonstrate that the 
embedded Co atoms form four covalent Co-Ge bonds, resulting in a Co 4+ valence state and a 3d 5 
electronic configuration. Calculated STM images are in perfect agreement with the experimental 
atomic resolution STM images for the broad range of applied tunneling voltages. 



PACS numbers: 68.43.Fg, 68.47.Fg, 68.37.Ef, 81.07.Ta 



I. INTRODUCTION 

The continuous miniaturization of electronic circuits 
has resulted in the emergence of novel classes of nanome- 
ter size devices that rel y on the quantum-mechanical na- 
ture of charge carriers!^ Examples of state-of-the-art 
nanodevices can be found in molecular electronics 3 ^ and 
spintronicsPThe dimensions of the active elements, con- 
nections and separations are now being reduced to the 
order of a few atomic rows and, in the ultimate limit, de- 
vices may be built up using atomic size elements^ that 
are connected by atomic nanowires. However, to fulfill 
the demands related to the never ceasing development of 
electronics, novel materials with electronic properties su- 
perior to that of the currently used silicon are required. 
Among all candidates germanium is consi dere d as one of 
the most promising alternative material s] 9 * 10 * because it 
allows higher switching speeds due to a lower effective 
hole mass and a higher electron and hole drift mobil- 
ityfCQ This makes germanium ideally suited for use in ul- 
trafast complementary metal-oxide-semiconductor tech- 
nology, in particular for metal-oxide-semiconductor field- 
effect transistors^"^ and band-to-band tunneling field- 
effect transistors.^ For this purpose detailed investiga- 
tions of the electronic properties of dopants/defects or 
metal alloys in Ge crystals are obviously required.^ 

One of the major challenges for future nanoelec- 
tronic applications is the controlled preparation of low- 
dimensional structures on semiconductor surfaces, e.g. 
quantum dots 1 and quantum wires ! 18 " 20 ! Due to their 
broad range of electronic and magnetic properties, such 
nanostructures are ideal model systems for the funda- 
mental study of low-dimensional physics as well as for 
the exploration of new device concepts^ th at also ex- 
ploit the spin character of the charge carriers!- 6 * 22 * 23 ! 



Within this context, deposition of metal atoms on 
Ge surfaces has attracted considerable scientific inter- 
est during recent years, since it was found that atoms 
of different materials self-organize into different types 
of nanostructures after deposition on Gc. It has been 
demonstrated that deposited Mn atoms do not coag- 
ulate on Ge(lll)c2x8 surfaces in the initial adsorp- 
tion sta ge, yi elding zero-dimensional (OD) str ucture s on 
Ge(lll)PMl On the other hand, PtBSBffl ArP 32 ! and 
Sn atoms' spontaneously form one- dimen sional (ID) 
atomic chains on Ge(001), whereas P d 34 ' 35 ' and Ag^ED 
atoms favor the formation of two-dimensional (2D) or 
three-dimensional (3D) particles on Ge(OOl). 

The emerging field of spintronics requires 
(self) assembly of nanostructures with well-defined 
magnetic properties on semiconducting surfaces. Due 
to the high spin polarization of the charge carriers 
near the Fermi level, Co is one of the most important 

elements used in magnetic recording media as well 

as in giant magnetoresistance devicesPO Recently, the 
electronic and magnetic behavior of ultrathin (< 5 
monolayers) Co/G<PHM] and Co/Ag/GePM] n i ms has 
been investigated. The initial adsorption stage of single 
Co atoms on Gc surfaces has not been studied so 
far. Thorough knowledge of the formation process of 
the Co/Ge interface during the first adsorption stages 
is, however, of crucial technological and fundamental 
interest. 

Here, we present a comprehensive study of the initial 
growth stages of Co on 2x1 reconstructed Ge(lll) by 
means of low-temperature (LT) scanning tunneling mi- 
croscopy (STM) and spectroscopy (STS), combined with 
first-principles density functional theory (DFT) calcula- 
tions within the local-density approximation. STM and 
STS are ideal tools to investigate with high spatial and 



energy resolution the surface reconstruction and the local 
electronic properties of the Ge(lll) surface after adsorp- 
tion of individual Co atoms. DFT calculations on the 
other hand allow us to predict the electronic properties 
of systems of up to thousands of atoms in size. High res- 
olution STM/STS combined with DFT hence provides 
a powerful tool for the investigation of atomic size sys- 
tems. Here, we report on the first experimental observa- 
tion of "noninvasive embedding" (i.e. without destroy- 
ing the surface reconstruction) of individual Co atoms in 
the Ge(lll)2xl surface and on the formation of larger 
Co/Ge intermixed layers after Co deposition on the cold 
(^sample < 80 K) Ge(lll)2xl surface. The location of an 
individual Co atom in the Ge surface, its influence on the 
surrounding Ge atoms and the resulting electronic prop- 
erties are systematically investigated. Voltage-dependent 
STM imaging reveals a highly anisotropic electronic per- 
turbation of the Ge surface surrounding the Co atom, 
which is accompanied by pronounced ID confinement 
along the 7r-bonded chains. Our experimental findings 
are well explained by the detailed DFT calculations. 



■*■ n-bonded chains 




[011] 



[211] 



[111] 



[211] 



FIG. 1. (Color online) (a) Schematic side view of the chain- 
left isomer of the 2x1 reconstruction of the Ge(lll) surface 
according to the original Pandey 7r-bonded chain model (dot- 
ted lines) and including the effect of buckling (solid lines) .^ 
The 7-member and 5-member Ge rings of the surface recon- 
struction are indicated by the numbers 7 and 5, respectively, 
(b) Schematic top view of "zigzag" chain structure of the three 
top layers in (a). The dashed frame indicates the surface unit 
cell, while the arrows with dotted lines indicate the "zigzag" 
structure of the upper 7r-bonded chain along [Oil] direction. 



II. INSTRUMENTATION 

STM and STS measurements were performed with a 
LT STM setup (Omicron Nanotechnology) , operating at 
a base pressure in the 10 -11 mbar range. All data are 
acquired at T samp i c ~ 4.5 K. Electrochemically etched W 
tips were cleaned in situ by repeated flashing well above 
1800 K to remove the surface oxide layer and any ad- 
ditional contamination. The tip quality was routinely 
checked by acquiring atomic resolution images o f the 
"herringbone" reconstruction of the Au(lll) surfacePZHSl 
STM topographic imaging was performed in constant 
current mode. The tunneling voltages Vt indicated in the 
text and figure captions are with respect to the sample 
(the STM tip is virtually grounded). Image processing 
was performed by Nanotec WSxMPSI 

Ge single crystals with a resistivity of pbuik — 0.2 51cm 
are doped with Ga at a doping level of ric a = 1 to 
2 x 10 16 cm~ 3 , resulting in p-type bulk conductivity. 
4 x 1.5 x 0.8 mm 3 Ge bars, with their long axis aligned 
with the (111) direction, were cleaved in situ at room 
temperature in the sample preparation chamber at a 
pressure of around 5 x 10 -11 mbar. The freshly cleaved 
samples were transferred within about 5 minutes to the 
STM measurement chamber. The pressure in the STM 
measurement chamber was about 4 x 10~ 12 mbar during 
the LT STM measurements. Under these conditions the 
cleaved Ge surfaces were observed to retain their clean- 
liness for 5 to 7 days. This way, we have investigated 7 
freshly cleaved Ge(lll) crystals. 

After checking the freshly cleaved Ge(lll)2xl sample 
in the STM measurement chamber, 0.02 to 0.04 mono- 
layers of Co are deposited on the cold Ge(l 11)2x1 sur- 
face (T samp i c < 80 K) in the sample preparation chamber. 
Deposition was achieved by evaporation from a high pu- 



rity Co (99.9996%) rod with an e-beam evaporator, at 
pressures below 10~ 10 mbar and at a low deposition rate 
of around 0.007 ± 0.001 monolayers (MLs) per seconds. 
After Co deposition, the Ge(lll)2xl sample is trans- 
ferred immediately to the LT STM measurement cham- 
ber. Overall transport time, including deposition time, 
was around 25 minutes. Here we focus on the results ob- 
tained on 4 different Co/Ge(lll)2xl samples for which 
2.1 /im 2 atomic resolution STM topography images were 
recorded and analyzed. 



III. EXPERIMENTAL RESULTS 

A. Topography and electronic structure of the 
freshly cleaved Ge(lll)2xl surface 

The 2x1 reconstruction of the Ge(lll) surface con- 
sists of 7r-bonded chains of Ge atoms running in the [011] 
direction.E2HSH Only every other (upper chain) row can 
be imaged by STM.^S The surface unit cell contains two 
atoms, both having one dangling bond. This dangling 
bond is responsible for 7r-bonding along the upper surface 
chain rows (see Fig. [I]) . In the original Pandey geome- 
trypU the two upper atoms of the 7-member ring have 
the same height and form "zigzag" chains along the [011] 
direction. However, due to buckling, one of these two 
atoms (up-atom) is shifted upwards (out of the surface) 
while the other (down-atom) is shifted downwards (into 
the surface) , as illustrated in Fig. [T] The occupied sur- 
face states are mainly localized on the up-atom, while the 
empty surface states are mainly localized on the down- 
atom. Consequently, the bonding surface states band 
7Tvb derived from the up-atom orbital is filled, while the 
anti-bonding surface states band 7r£. B derived from the 
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FIG. 2. (Color online) (a)-(c) Typical large scale STM to- 
pography images of the freshly cleaved Ge(lll)2xl surface 
(Vt = +1.0 V, It = 35 pA). High resolution (d) empty and 
(e) filled states STM images of the same area, recorded at the 
indicated tunneling voltage Vt and at 7 t = 0.9 nA and 3.0 nA, 
respectively. 



down-atom orbital is empty. 

In Fig. [2] we present typical large scale [(a)-(c)] and 
high resolution [(d)-(e)] STM topography images of the 
clean Ge(lll)2xl surface. Large atomically flat terraces 
up to 10 5 nm 2 can be observed, which are separated from 
each other by monatomic steps (MASs). It can be ob- 
served in Fig. [2] (b) that the Ge(lll)2xl surface con- 
sists of different types of domains with slightly different 
atomic arrangement. 54 This is related to the threefold 
rotational symmetry of the surface. The domains are 
found to be separated by two different types of domain 
boundaries (DBs) [see Fig. [2] (b)]. At the first type of 
DB, referred to as a type-A DB following the terminol- 
ogy used in Ref. |54j . the atomic rows at the opposite 
sides of the DB are rotated by an angle n/3. The second 
type DB, the so-called anti-phase DB or typc-B DB, [54] 
is formed due to a shift of the 7r-bonded chain rows in the 
[211] direction by half a unit cell. We found that most 
DBs are of type-B and that the type-A DBs often exhibit 
(local) disorder! 5 - 5 -' 

In addition, two types of MASs can be observed [see 
Figs. [2] (a) and (c)]. With respect to the upper terrace 
one type of MASs, hereafter referred to as type-A, is 
oblique to the 7r-bonded chain rows on the terrace. The 
second type of MASs, hereafter referred to as type-B, is 
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FIG. 3. Normalized conductance spectrum recorded on a de- 
fect free area of the freshly cleaved Ge(lll)2xl surface. The 
main energy bands are indicated by gray rectangles at the bot- 
tom of the hgure (see text for more details). Inset: Current- 
voltage (I-V) characteristic close to the Fermi energy Ef. 



parallel to the 7r-bonded chain rows. One should note 
that Fig. [2] (a) exhibits pronounced Moire fringes that 
run along one direction and become visible because of 
the large size of the imagepS while the 7r-bonded chain 
rows are not visible on the STM image and their direction 
is specified by arrows for each of the terraces. We find 
that Ge adatoms are often present at the MASs, both on 
the upper and on the lower terrace, except on the upper 
terrace of type-A MASs. At these terraces the surface is 
either locally distorted or a 2x4 or c2x8 surface recon- 
struction occurs [Fig. (c)].^ The Ge surface adatoms 
are probably created upon cleavage at room temperature, 
after which the adatoms can migrate along 7r-bonded 
chain rows to the MAS regions. Furthermore, (individ- 
ual) Ge adatoms can be frequently observed on atom- 
ically flat Ge(lll)2xl terraces as well above a charged 
subsurface Ga impurity [see label (1) in Fig.|2](c)]. These 
adatoms are well separated from each other and their 
number is in good agreement with the low doping level 
of our Ge samples. 

In Fig. [3] we present a typical normalized conductance 
spectrum recorded at the p-type Ge(lll)2xl surface. 
The main energy bands are indicated by gray rectan- 
gles at the bottom of Fig. [3] In the spectrum the large 
peak around 0.19 eV can be assigned to the onset of the 
unoccupied surface states conduction band itq B related 
to the uppers-bonded chains of the Ge(lll)2xl surface 
(see Fig. n])! 57 * 58 ! Two energy gaps can be discerned: A 
narrow gap of about 0.19 eV and a wide gap of about 
0.74 eV. The latter corresponds to the forbidden energy 
gap of the projected bulk band structure of the Ge(lll) 
surface at low temperature. On the other hand, the 
narrow gap corresponds to the energy gap between the 



filled bulk valence band (VB) and the unoccupied surface 
states conduction band (7T ^ R ). Using STS^l as we n as 



photoemission experimentspSHEU the surface states band 
gap has been determined before, yielding a gap value of 
0.54 ± 0.04 eV. From Fig. [3] it can be concluded that the 
high resolution STM images in Fig. [2] (e) and Fig. [2] (d) 
were obtained at tunneling voltages near the top of VB 
and the bottom of tTq B , respectively. 



B. Adsorption of Co atoms on Ge(lll)2xl 

In Figs. El (a) and (b) we present two typical large 
scale STM topography images of the Ge(lll)2xl surface 
after Co atom deposition corresponding to a coverage of 
0.032 ML. Three different kinds of structures are formed 
after Co deposition: (i) Co/Ge intermixing layers (ILs) 
[indicated by the two arrows with label 1 in Fig. 4] (a)], 
(ii) Co clusters consisting of multiple Co atoms [indicated 
by the two arrows with label 2 in Fig. [4] (b)] and (iii) 
well separated individual Co atoms [indicated by the two 
arrows with label 3 in Fig. [41(b)]. 

The Co/Ge ILs are formed due to the consecutive ac- 
cumulation of Co atoms at surface/subsurface defects, 
e.g. DBs and MASs. 62 As indicated by the results of our 
DFT calculations that are presented below, Co atoms are 
able to migrate along the 7r-bonded chain rows, despite 
the low temperature of the sample during Co deposition 
(T s < 80 K). Co/Ge ILs are found both on the upper and 
lower terraces at type-A MASs. At the type-B MASs, 
formation of a Co/Ge IL occurs only on the lower ter- 
race. Near DBs Co/Ge ILs are observed on both sides of 
the type-A and type-B DBs. Co/Ge ILs can be found on 
atomically flat terraces as well, far away from any DBs 
and MASs. Since Co free atomic size defects, including in 
particular Ga subsurface impurities,^ can no longer be 
observed after Co deposition, this suggest that these de- 
fects act as nucleation centers for the formation of Co/Ge 
ILs. The amount of Co/Ge ILs formed on atomically flat 
terraces roughly scales with the amount of atomic scale 
defects that is observed prior Co deposition. A more de- 
tailed discussion on the formation of Co/Ge ILs will be 
presented elsewhere!^! 

Only a small fraction of the Co atoms coagulates into 
small Co clusters. A larger fraction of the Co atoms 
remains under the form of individual atoms after deposi- 
tion. At tunneling voltages above 0.7 V these individual 
atoms appear as bright protrusions located on the upper 
7r-bonded chain rows. At other tunneling voltages the 
Co atoms are observed differently. This voltage depen- 
dence is discussed in more detail in Section fill CI below. 
High resolution STM topography images of individual Co 
atoms are presented in Figs.E (c) and (d). The amount of 
individual Co atoms that can be inferred from Figs. El (a) 
and (b) is 0.005 ± 0.002 ML. 13 ± 5% of the deposited 
amount of Co atoms is observed as individual Co atoms, 
while 87 ± 5% contributes to the formation of Co/Ge ILs 
and Co clusters. It is important to already note here that 



# 




FIG. 4. (Color online) (a), (b) STM topography images of the 
Ge(lll)2xl surface after deposition of 0.032 ML of Co (V t = 
+ 1.0 V, It = 15 pA). (c), (d) High resolution STM images of 
the areas confined by the dashed rectangle and by the dotted 
rectangle in (b), respectively (Vt = +0.9 V, It — 10 pA for (c) 
and V t = +0.80 V, h = 300 pA for (d)). 



the individual Co atoms are actually not on top of the Ge 
surface, as will be demonstrated in detail in Section IV 
below by comparing the experimental STM images to 
simulated STM images based on DFT calculations. Ac- 
cording to the DFT calculations, individual Co atoms 
penetrate into the Ge surface and reside in between the 
3 rd and the 4 th atomic layer (AL), in a quasi-stable posi- 
tion (at low temperatures) inside the 7-member Ge ring 
of the 2x1 reconstruction (see Fig. [I]). This "embed- 



ding" is found to influence the local electronic structure, 
but does not give rise to a modified surface reconstruc- 
tion, as can be seen in Figs. El (c) and (d). We therefore 
refer to this embedding as "noninvasive" . 

The low temperature of the Ge substrate during Co 
deposition appears to be crucial for obtaining individ- 
ual, well separated Co atoms. As mentioned above, in 
spite of the low substrate temperature, the larger fraction 
(around 87%) of the Co atoms still exhibits sufficient sur- 
face mobility to migrate to defects where Co/Ge ILs are 
formed [see Fig. [41(a)]. The residual fraction of deposited 
Co atoms (around 13%) remains confined to defect free 
atomically flat terraces of the 2x1 surface. Recently, we 
reported that these Co atoms diffuse from their quasi- 
stable sites to surface/subsurface defects as well after 
warming up the sample to room temperature.^ 2 How- 
ever, as long as the sample remains at low temperatures, 
the embedded individual Co atoms remain immobile dur- 
ing the experiments in the investigated —1.5 to + 1.5 V 
voltage range. 

As illustrated in Fig. EJ (d) , the "zigzag" structure of 
the 2x1 reconstruction of the low doped p-type Ge(lll) 
surface (see Fig. IT]) becomes observable for a limited 
range of tunneling voltages V z [ gzag — 0.85 ± 0.07 V. 
From Fig. [3] we can conclude that both the unoccupied 
surface states itq B (wave functions that are mainly local- 
ized on the down-atoms) and the unoccupied states at the 
bottom of conduction band (CB) (wave functions which 
are partially localized on the wp-atoma^) become avail- 
able for tunneling within this voltage range. Although 
the wave functions on the Mp-atoms have a smaller am- 
plitude, the higher position of these atoms implies that 
they appear more prominently in the STM images than 
the down-atoms when the applied tunneling voltage in- 
creases. Within the V z - lgzag voltage range (this range is 
marked by the gray dashed area in Fig. pi), the contribu- 
tion of ttq B to the tunneling current remains nearly con- 
stant, while the contribution from CB rapidly increases 
with increasing Vt up to 1.0 V. Hence, a balance exists 
between tunneling into up- and down-atoms within the 
^zigzag voltage range, which implies that both the up- 
and down-atoms of the 7r-bonded chain rows are visu- 
alized in constant-current STM images [see Fig. EJ (d)]. 
We found that the precise value of V z i gzag depends on 
the doping level as well as on the semiconductor type. 
E.g., the 2x1 reconstruction of the (111) surface of 
heavily doped n-type Ge(lll) (phosphorus doping level 
rip = 1 x 10 19 cm~ 3 ) reveals zigzag chains around 0.55 V. 
Previously, it has been reported by Trapmann et al. that 
zigzag chains appear around 0.8 V for n-type Si(lll)2xl 
surfaces P^l The precise value of V zlgzag is therefore a char- 
acteristic feature of the semiconductor surface under in- 
vestigation. 

From Figs. EJ (a)-(d) it becomes clear that all individ- 
ual Co atoms occupy identical positions at the Ge surface: 
Their appearance is symmetric with respect to the [211] 
direction [the symmetry axis is drawn as a black dash- 
dotted line in Fig. Pfl (d)] , while they appear asymmetric 
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FIG. 5. (Color online) High resolution STM topography 
image of 5 Co atoms embedded in the Ge(l 11)2x1 surface 
(K = +0.8V, I t =200pA). 



with respect to the [011] direction. The latter can be 
related to the asymmetry of the 2x1 surface reconstruc- 
tion (due to buckling) along this direction. The relevant 
crystallographic directions are indicated in Fig. EJ (c) and 
are the same for Figs.H(b) and (d). The symmetry and 
asymmetry of the Co atoms with respect to the [211] 
and [Oil] directions, respectively, becomes most clearly 
resolved for tunneling voltages near V z \ gzag [see Fig. EJ (d) 
and Fig. pi. The bright protrusion related to an indi- 
vidual Co atom occupies about 2 surface unit cells along 
the 7r-bonded chain rows. Moreover, in Fig. EJ (d), dis- 
turbance of the zigzag atomic structure along the [Oil] 
direction can be observed near the Co atom over a dis- 
tance < 3 unit cell periods. In Fig. EJ (d) the edge of the 
unperturbed zigzag 7r-bonded chain on the right hand 
side of the Co atom is marked by a long dotted white 
line, forming an angle a s=s 35° with the [211] sym- 
metry axis. Closer to the Co atom related protrusion, 
a brighter/higher feature [marked by the short dotted 
white line and indicated by the label 2 in Fig. EJ (d)] is 
visible, also making an angle close to a with the [211] 
symmetry axis. The same disturbance is observed to the 
left of the Co atom (mirror symmetry with respect to 
the [211] direction). Identical characteristic features are 
observed for all investigated individual Co atoms on dif- 
ferently oriented Ge(lll)2xl domains. 

In addition to the disturbance of the zigzag structure, 
an asymmetry of the electronic structure of the neigh- 
boring 7r-bonded chain rows (with respect to the [Oil] 
direction) near an embedded Co atom is observed. This 
asymmetry is most clearly resolved in maps of the local 
density of states (LDOS) (not shown here). It can be ob- 
served as well in STM topography images around V z - lgzag : 
The lower 7r-bonded chain row in Fig. RJ (d) (indicated by 
the green solid arrow) appears brighter when compared 
to subsequent rows, while the upper 7r-bonded chain row 




FIG. 6. (Color online) (al)-(a6) Empty states and (bl)-(b6) filled states STM topography images of 3 single Co atoms embedded 
in the Ge(lll)2xl surface. Image sizes are 15.0 x 11.2nm 2 . The tunneling voltage Vt is indicated for each image. Tunneling 
current I t is 180pA, 300pA, 320pA, 340pA, 340pA, and 340pA for (al,bl) to (a6,b6), respectively. 



(indicated by the green dashed arrow) does not appear to 
be influenced by the embedded Co atom. This asymmet- 
ric perturbation and the disturbance of the zigzag struc- 
ture of the 7r-bonded chain rows always occur simultane- 
ously and were observed with various STM tips. Figure[5] 
presents an STM topography image of 5 embedded single 
Co atoms that all induce a similar difference between the 
Co neighbor chains. As we will demonstrate below when 
discussing the results of our DFT calculations, this addi- 
tional asymmetry of the Ge (111)2x1 surface induced by 
the Co atoms allows us to determine the crystallographic 
[211] direction of the different Ge(lll)2xl domains, lead- 
ing to the conclusion that all investigated Ge(lll)2xl 
surfaces consist of domains with 7r-bonded chain-left iso- 
mers (negative buckling) exclusively. 

There appears no significant electronic interaction be- 
tween two embedded single Co atoms that reside in the 
same 7r-bonded chain row in close vicinity (down to a dis- 
tance of 3 unit cells along the [Oil] direction) in the inves- 
tigated voltage range from —1.5 to +1.5 V [see Figs. [4] (a) 
and (b), as well as Co atoms (3) and (4) in Fig. [5J. On 
the other hand, an electronic interaction effect can be ob- 
served for Co atoms that reside in neighboring 7r-bonded 
chain rows at a distance (in the 7t-bonded chain row di- 
rection) smaller than ±3 unit cells. This is illustrated 
in Fig. p)J Indeed, the lower Co atom (3), which is lo- 
cated on a 7r-bonded chain row that is perturbed by the 
upper Co atom (2), exhibits modified electronic proper- 
ties. More precisely, Co atom (3) appears more bright in 
Fig. [51 while Co atom (2) appears similar to the other 
Co atoms. Note that the modified electronic properties 
of Co atom (3) are caused by Co atom (2) only and are 
not related to presence of the neighboring Co atom (4) 



that is located in the same 7r-bonded chain row. 



C. Voltage dependent STM investigation of single 
Co atoms embedded in Ge(lll)2xl 

In this section we focus on the novel Co induced elec- 
tronic features by careful comparison to the electronic 
properties of the clean Ge(lll)2xl surface (see Sec- 
tion 



III A). First of all it must be noted that deposition 



of (a small amount of) Co atoms does not change the 
electronic properties of the defect free Ge(lll)2xl sur- 
face: The characteristic peaks in the STS spectra of the 
clean Ge(lll) surface prior to Co deposition (Fig. [3]) are 
still observed at the 2x1 reconstructed surface after Co 
deposition. Figure k$\ presents a series of (a) empty and 
(b) filled states STM topography images of 3 individual, 
well separated Co atoms embedded in the reconstructed 
Ge(lll)2xl surface. Images are recorded at the same 
location and with the same tip for a broad range of tun- 
neling voltages V t between 1.4 V and —1.5 V. The series 
of images reveals a pronounced voltage dependence for 
both the clean p-Ge(lll)2xl surface and the electronic 
influence of the embedded Co atoms on the 7r-bonded 
chain rows. 

In the empty state regime, at high voltages above 
1.0V [Figs. [6] (a5) and (a6)], topography is dominated 
by the reconstruction lines of the 7r-bonded chain rows 
in the [Oil] direction. Individual Co atoms appear 
as bright protrusions, located directly on the upper tt- 
bonded chain row and extending over 2 unit cells of the 
2x1 reconstruction. Maximum contrast of the Co related 
features is observed around V t = 0.93 V and above this 
voltage the contrast again decreases. This is because the 



contribution of the ttq B surface states to the tunneling 
current remains approximately constant, while the con- 
tribution of the bulk CB states increases with increasing 
tunneling voltage (see Fig. [3]). At lower voltages around 
Kdgzag = 0.80 V [Fig. [6] (a3j] the zigzag atomic structure 
discussed above emerges. Around V t = 0.60 V contrast 
of the Ge atomic corrugation and the Co related protru- 
sions become similar and Co atoms can mainly be dis- 
cerned by the locally induced perturbation of the 2x1 
surface reconstruction [Fig. [6] (a2)]. In the empty states 
regime below 0.60 V [Figs. pRal) and (a2)] extra corru- 
gation appears in the STM topography along the [211] 
direction. Close to the Fermi level Ep the strength of 
this extra corrugation along the [211] direction becomes 
comparable to the corrugation along the [011] direction 
[Fig. [6] (al), also see Fig. [2] (d)]. Co atoms appear as 
centro-symmetric striped depressions along the 7r-bonded 
chain rows. These depressions exhibit a local minimum 
directly above the embedded Co atom and gradually fade 
away with increasing distance (up to 4nm) from the local 
minimum [Fig. ro] (al)]. 

In the filled states regime, at voltages close to Ep, to- 
pography is again dominated by corrugation of the n- 
bonded chain rows [Fig.[6](bl), also see Fig. [21(e)]. Max- 
ima of the atomic corrugation are related to the highest 
filled bulk VB states that are localized on the Ge up- 
atoms (see Fig. [I]). Co atoms appear as bright protru- 
sions located on the 7r-bonded chain rows [Fig. [6] (bl)]. 
Below —0.7V extra corrugation emerges along the [211] 
direction [Fig. p^ (b4)] and persists down to around 
— 1.2 V. Below —1.2 V topography becomes again com- 
pletely dominated by the 7r-bonded chain rows along the 
[011] direction [Fig. p] (b6)]. Here, Co atoms appear as 
atomic size vacancies in the upper 7r-bonded chain rows 
of the 2x1 surface reconstruction [Figs. p] (b5) and (b6)]. 

It is clear from Figs.[6](al) and (b2)-(b5) that the em- 
bedded Co atom induces a localized ID perturbation of 
the LDOS along the 7r-bonded chain rows. STM images 
recorded with a tunneling voltage near the edge of the 
surface states bands reveal highly anisotropic scattering 
of electrons and screening effects with ID confinement to 
the 7r-bonded chains. ID perturbations have previously 
been reported for Si(lll)2x]PESI and Ge(lll)2xl^l 
surfaces. Here, the perturbation near the atomic size 
defect has a pronounced ID shaped and extends up to 
6 nm along the 7r-bonded chain rows, while the width of 
the perturbation remains limited to one period (0.69 nm) 
of the 2x1 reconstruction in the [211] direction. In the 
empty states regime, the ID screening effects appear as 
depressions and become most pronounced near 0.23 eV, 
corresponding to the bottom of the surface states band 
ttq B [Fig.[6](al)]. In the filled states regime, the ID elec- 
tron scattering effects are observed as protrusions and 
are most clearly seen near —0.66 V, corresponding to the 
top of the surface states band ttvb [Fig.|6](b3)]. In both 
regimes the ID perturbations exhibit identical mirror-like 
symmetry with respect to the [211] direction, similar to 
our recent observations for atomic size surface impurities 



on Ge(lll)P 

In summary, we can state that individual Co atoms 
embedded in Ge(lll)2xl surfaces exhibit the following 
general properties: 

• Individual Co atoms penetrate into the cold (T s < 
80 K) Ge(lll)2xl surface after deposition (see Sec- 
tion 



IV below for more details); 



• the 2x1 reconstruction is preserved after the non- 
invasive embedding of a Co atom; 

• the embedded Co atoms occupy identical positions 
in the Ge(lll) surface and exhibit an identical volt- 
age dependence of the STM topography images; 

• the embedded Co atoms exhibit a clear symmetry 
with respect to the [211] direction, while they ex- 
hibit a clear asymmetry with respect to the [Oil] 
direction; 

• the embedded Co atoms induce highly anisotropic 
scattering of electrons, which is accompanied by 
screening effects with ID confinement along the tt- 
bonded chain rows. 



IV. DISCUSSION AND COMPARISON WITH 
DFT CALCULATIONS 

A. DFT model of the Ge(lll)2xl surface 

The 2x1 reconstruction of cleaved Si and Ge (111) 
surfaces^lSIHZni j s we \\ described by the commonly used 
Pandcy 7r-bonded chain model described abovepD The 
Pandey chain geometry leads to a strong coupling of 
the dangling-bond orbitals along the chain, while the 
coupling between the chains is much weaker. This ge- 
ometry does not take into account the effects of buck- 
ling, so that the two uppermost surface Ge atoms in the 
7- member rings are at the same height!^ Northrup et 
al. predicted buckling of the two uppermost atoms by 
about 0.8A (see Fig. [l]).^2l This buckling further reduces 
the surface energy of the system, yielding two different 
isomers commonly referred to as 7r-bonded "chain-left" 
and "chain-right" isomers^ (see Fig. Uj . Relying on first 
principle calculations within the computational accuracy 
that could be achieved at that time (1991), it was found 
that the chain-left isomer is 6 meV per lxl surface cell 
lower in energy when compared to the chain-right isomer. 
It has be en confirmed theoretically^ and experimen- 
tall} * 63 ! 73 -! that the chain-left isomer is indeed the dom- 
inant isomer at the Ge(lll)2xl surface. However, be- 
cause of the very small energy difference between the two 
different isomers, coexistence of both surface reconstruc- 
tions and hence the possible existence of Ge (111)2x1 
multi-domain surfaces cannot be totally excluded. 

Our theoretical investigation of the noninvasive em- 
bedding of Co atoms in the Ge(lll)2xl surface was 
performed based on DFT^l within the local density ap- 
proximation (LD A). Calculations were performed with 
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FIG. 7. (Color online) The final position of the 9 topmost Ge 
layers of the Ge(lll)2xl surface reconstruction: (a) Chain- 
left and (b) chain-right isomers of the Pandey 7r-bonded chain 
model with buckling (side view), (c) Top 3 surface atomic 
layers of the 4 x 2 SC chain-left Ge (111)2x1 isomer model 
(top view), (d) Surface Brillouin zone and relevant directions 
of the Ge(lll)2xl surface. 



the SIESTA packagep^S which relies on the expan- 
sion of the Kohn-Sham orbitals by linear combination 
of pseudo-atomic orbitals. In all calculations a doublc- 
zeta basis set with polarization was used. The core elec- 
trons were implicitly treated by using norm-conserving 
Trouiller-Martins pseudopotentials^ with the following 
electronic configuration of the elements: H Is 1 , Ge (Ar 
3d 10 ) 4s 2 4p 2 and Co (Ar) 4s 2 3d 7 , where the core con- 
figurations are indicated between parentheses. A cutoff 
energy of 200 Ry was introduced for the grid integration, 
ensuring convergence of the total energy of the system 
within typically 0.1 meV. 

Our calculations are performed in three stages. First, 
we determine the surface equilibrium geometry for the 
chain-left and chain-right isomers within one surface unit 
cell (SUC) using conjugate gradient (CG) geometry op- 
timization. Second, based on the Ge(lll)2xl SUC, a 
new n x m supercell (SC) with one Co atom above the 
surface is constructed and the equilibrium configurations 
of Co/Ge(lll)2x 1 are determined relying on CG geome- 
try optimization. Third, the electronic structure of a still 
larger SC is calculated, which is then used to perform the 



DFT based STM topography image simulations. 

In order to model the Ge(lll)2xl surface, a reduced 
square SUC (unreconstructed) with (ao\/2, arj\/6) lattice 
vectors is used (ao = 2.8205A is half the optimized bulk 
Ge lattice constant). This SUC consists of a slab of 26 
Ge atomic layers, of which one atomic layer is saturated 
by hydrogen atoms (52 Ge atoms and 2 H atoms per 
SUC). In the CG geometry optimization the 14 topmost 
atomic Ge layers are allowed to move, while the 12 layers 
of Ge and H atoms are frozen at the ideal (bulk) posi- 
tions. Using 130 fc-points within the surface Brillouin 
zone, the atoms are relaxed until all atomic forces acting 
on the released atoms are smaller than 3 meV/ A, and the 
remaining numerical error in the total energy is smaller 
than 0.1 meV for each optimization step. 

The final positions of the 9 topmost Ge layers yield- 
ing the Ge(lll)2xl surface reconstruction are presented 
in Figs. [7] (a) and (c) (chain-left isomers) and Fig. [7] (b) 
(chain- right isomers) . The coordinates of the Ge atoms 
below the 7 th atomic layer are found to change only 
slightly during the CG geometry optimization. For the 
considered SUC the two different Ge isomers are both 
possible for the formation of the 2x1 reconstruction. De- 
pending on the initial conditions, either the 5-member or 
the 7-member Ge ring of the reconstruction is present in 
the 2x 1 SUC. This is illustrated in Figs. [7] (a) and (c) for 
the chain-left isomer. The black solid bars in Figs. [7] (a) 
and (c) comprise the 5-member ring SUC, while the red 
dashed bars comprise the 7-member ring SUC. The same 
can be done for the chain-right isomer, as illustrated in 
Fig. u\ (b). For the 5-member ring the upper 7r-bonded 
chain rows are formed at the joint of the SUC in the [211] 
direction and are referred to as "chain-left /right (1)" 
hereafter [see, e.g., the solid colored atoms in Fig[7](b)]. 
For the 7-member ring, the upper 7r-bonded chain rows 
are formed inside the SUC and are referred to as "chain- 
left/right (2)" hereafter [see, e.g., the solid colored atoms 
in Fig[7](a)]. We verified that the chain-left /right (1) and 
chain-left/right (2) selection exhibit identical electronic 
properties for periodical boundary conditions. Depend- 
ing on the position of the adsorbed Co atom with respect 
to the upper 7r-bonded chain rows in the n x m SC, se- 
lection (1) or (2) was chosen. Buckling distances of the 
chain-left and chain-right isomers are 0.83A, and - 0.80A, 
respectively. The chain-right isomer has a total surface 
energy that is 14meV/(2 x 1 SUC) higher than the total 
surface energy of the chain-left isomer, implying that the 
chain-left isomer should be the donrinant isomer for the 
Ge(lll)2xl surfacePEUm 



B. Co adsorption sites and energy decomposition 

To model the adsorption of a Co atom on the 
Ge(lll)2xl surface, we used the 9 topmost relaxed Ge 
layers for the chain-left and the chain-right isomer that 
were obtained following the procedure described above 
(also see Fig. [7| . 6 of these 9 atomic layers are allowed 
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to relax. The bottom side of the 3 fixed atomic layers 
is saturated by hydrogen atoms (16 Ge atoms and 2 H 
atoms per SUC; 16 A slab vacuum space separation). 

Geometry optimization was carried out for an enlarged 
4x2 SC for both the chain- left and the chain-right iso- 
mers. In Fig. [7] (c) we show the three topmost layers 
of the 4x2 SC (size is 15.95 A x 13.81 A, consisting of 
145 atoms) for the chain- left (1) isomer geometry. A Co 
atom was then located in front of the Ge (111)2x1 sur- 
face. Next, the relaxation of the Co atom was calculated 
by means of CG geometry optimization. CG geometry 
optimization for the 4x2 SC was carried out until all 
atomic forces acting on the released Co and Ge atoms 
were below 5 meV/A and until the numerical error on 
the total energy was smaller than 10~ 4 eV per SC for 
each optimization step. Using a variety of starting co- 
ordinates for the Co atom, multiple Co/Ge quasi-stable 
geometries were tested and their total energies were com- 
pared. Both chain-left and chain-right isomer geometries 
were used. In both cases multiple quasi-stable Co atom 
sites (with respect to an atomic force tolerance 3 meV/A) 
were found, including sites on the Ge(lll)2xl surface as 
well as underneath the Ge(lll)2xl surface, i.e., inside 
the big 7-member Ge ring. The subsurface sites were 
identified by using a location inside the 7-member ring 
as the starting location for the CG geometry optimiza- 
tion. On the other hand, the surface sites were identified 
by using an initial Co location above the Ge(lll)2x 1 sur- 
face, at heights in the 3 to 4 A range and at various initial 
coordinates in the (x, j/)-plane. This way, 9 quasi-stable 
sites were identified for the chain-left isomer, of which 6 
are located on the Ge surface and 3 are located inside 
the 7-member ring underneath the Ge surface. For the 
chain-right isomer 5 surface and 3 subsurface quasi-stable 
sites were identified. 

In Table |l] we present for the chain-left isomer geome- 
try and for all possible quasi-stable Co atom sites (labeled 
(1) to (9) in the first column) an overview of the calcu- 
lated difference AE S i tC s in total energy with respect to 
the total energy for the Co atom located at site (7). For 



TABLE I. Calculated difference A_E s it es in total energy for 
the Co atom located at the different sites (on top (1-6) and 
underneath (7-9) the Ge(lll)2xl surface) and the Co atom 
located at site (7). The calculated energy differences are for 
the chain-left isomer geometry and are given for 4x2, 8x4 
and 14x4 SCs. 



Co site number 


4x2 SC 


8x4 SC 


14x4 SC 


(1) 


1.791 eV 


moves to site (2) 


- 


(2) 


0.376 eV 


0.724eV 


1.250 eV 


(3) 


0.583 eV 


0.730 eV 


0.841 eV 


(4) 


0.966 eV 


1.199 eV 


1.350 eV 


(5) 


2.280 eV 


2.470 eV 


- 


(6) 


0.701 eV 


1.054 eV 


1.280 eV 


(7) 


OeV 


OeV 


OeV 


(8) 


0.038 eV 


0.129 eV 


0.462 eV 


(9) 


0.512 eV 


0.595 eV 


0.693 eV 



a 4x2 SC size (second column) the minimum energy was 
found to occur for the Co atom located at site (7). Simi- 
lar calculations were performed for the chain-right isomer 
geometry and for a 4x2 SC (data not shown). Again, the 
minimum energy was found for the case of a Co atom lo- 
cated inside the 7-member ring, but with somewhat lower 
energy gain when compared to the chain-left geometry. 

In order to reduce the influence (related to the peri- 
odic boundary conditions) of the restructuring induced 
by the Co atom within a 4x2 SC, similar calculations 
were performed for a 8x4 SC (size is 31.90A x 27.63 A, 
consisting of 577 atoms). The obtained values for the 
difference AE'sitos in total energy for all possible quasi- 
stable Co atom sites after CG geometry optimization are 
listed in the third column of Table [J] Similar to the case 
of the 4x2 SC, the energy difference is again with respect 
to the total energy for the Co atom at site (7). When us- 
ing this larger SC, a transition of the Co atom was found 
from site (1) to site (2). Also, site (9) was found to be 
energetically more favorable than site (2) when compared 
to the 4x2 SC. This can be explained by the occurrence 
of longer-range surface relaxations inside the larger SC. 
Remarkably, for the 8x4 and 14x4 SCs all Co sites inside 
the 7-member ring have a lower energy (see Table fil) with 
respect to the other sites, which is different from the 4x2 
SC. 

Finally, calculated energy values for a 14x4 SC (size 
is 55.82 A x 27.63 A, consisting of 1009 atoms) are pre- 
sented for the chain-left isomer geometry in the fourth 
column of Table III Similar calculations were performed 
for the chain-right isomer geometry (data not shown) and 
yielded nearly identical results, with again a somewhat 
lower energy gain with respect to the other sites when 
compared to the chain-left geometry. 



C. Prime location of the Co atom into Ge(lll)2xl 
surface 



In Fig. [8] (a) we present a 3D ball-and-stick model view 
of the relaxed Ge(lll)2xl surface, obtained for the 8x4 
SC with chain-left isomer geometry. The Co atom is 
located at the site of minimum energy, i.e. site (7) (see 
third column in Table ffl. A 3D front view and a side view 
of the Co atom inside the 7-member ring are presented 
in Fig. [8] (b) and in Fig. [8] (c) , respectively. Relaxation of 
the surface Ge up-atoms and down-atoms upon Co atom 
incorporation can be clearly observed. The black solid 
envelope line in Fig. [8] (c) reflects the variation of the z- 
coordinate of the center of the Ge up-aioms (magnified 
by a factor 2 for clarity). The upward shift of the Ge 
up-atom located directly above the Co atom is +0.43A, 
while the downward shift of the neighboring Ge up-atoms 
is — O.IOA. The Co induced shift of the z-coordinate of 
the Ge up- and down-atoms extends as far as ±3 peri- 
ods along the 7r-bonded chain row of the 2x1 reconstruc- 
tion, in agreement with our experimental observations 
[see Fig. El (d) and the related discussion in Section III B 
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FIG. 8. (Color online) Calculated minimum energy position 
of a single Co atom (yellow) at the Ge(lll)2xl surface. The 
Co atom is in between the 3 rd and 4 th atomic layer underneath 
the Ge surface, (a), (b) 3D view of the chain-left Ge isomer. 
The Co atom is located inside the 7-member Ge ring (viewed 
along the [Oil] direction), (c) Side view of the Co-host chain- 
left isomer (viewed along the [211] direction). Height varia- 
tions of the itp-atoms in the 7r-bonded chain row are indicated 
by the black envelope curve (magnified by a factor 2 for clar- 

ity). 



The black solid arrows in Fig. [8] (b) (left image) indi- 
cate the bulk Ge atoms that experienced the most signif- 
icant shift of their positions upon embedding of the Co 
atom. In Fig. M (b) (right image) it can be seen that the 
Ge atoms to the left of the Co atom (indicated by the 
black dashed arrows) remain unperturbed. This asym- 
metry of the geometry (and hence of the local electronic 
properties) of the Co/Ge(lll)2xl system along the [Oil] 
direction is in agreement with our experimental observa- 
tions: The STM topography images in Fig. [J] (d) and 
Fig. [5] also exhibit an asymmetry with respect to the 
[Oil] direction around an embedded Co atom. Finally, 
we want to stress once more that the embedding of Co 
atoms does not give rise to a novel Ge surface reconstruc- 
tion. Instead, the Ge atoms surrounding the Co atoms 
experience only small changes of their positions, which 
is accompanied by changes of the local electronic prop- 
erties as well. Experimentally, we also found that the 
2x1 reconstruction is maintained upon Co embedding, 
as becomes clear in Figs. [6] (al), (a2), (bl) and (b2). 
Calculations for the other Co atom sites in Table U re- 



veal that these sites lead to more drastic changes and 
in some cases even local destruction of the Ge(lll)2xl 
reconstruction. 



D. Embedding of a Co atom into the Ge(lll)2xl 
surface 



In this section we will discuss the two most feasible 
routes for the incorporation/penetration of a Co atom, 
which is initially above the Ge(lll)2xl surface, to the 
site (7), i.e. inside the 7-member ring of the Ge(lll)2xl 
surface. Embedding of deposited atoms into subsurface 
layers has already b een d emonstrated before for a SP^l 
surface and for a G d 81 * 82 * surface. Ge atoms have been 
found both experimentally and theoretically to penetrate 
into the 4 th subsurface layer of Si(lOO)^- when deposited 
at a substrate temperature of about 500° C. Similarly, it 
has been found that Si atoms deposited on Ge(100)2xl 
are able to move below the Ge surface at room tempera- 
ture!^ Finally, the formation of Co/Ge intermixing lay- 
ers of up to 3 MLs thick after deposition of Co atoms on 
room temperature GeQll) substrates has been demon- 
strated experimentally!^! 

As a possible starting site for penetration into the 7- 
member ring, the energetically two most favorable quasi- 
stable Co atom sites on top of the Ge(lll)2xl surface 
are considered: Site (2) and site (3) (see Table [I]). A 3D 
visualization of both configurations is presented in Fig. [9] 
The crystallographic directions of Fig. [9] are identical to 
those of Fig. [8] (a) (except for a rotation of the viewpoint 
angle with respect to the [111] direction). For the sites (2) 
and (3) the Co atoms are located in between two neigh- 
boring upper 7r-bonded chain rows, at the center of the 
hexagonal 6-member Ge ring on the left [Fig. [9] (a)] and 
right hand side [Fig. [91(b)] of the central upper 7r-bonded 
chain row. A top view of these sites at the Ge(l 11)2x1 
surface is presented in Fig. [To] (b) . The Co atom at site 
(2) has one neighboring Ge down-atom and two neighbor- 
ing Ge wp-atoms, while the Co atom at site (3) has one 
neighboring Ge up-atom and two neighboring Ge down- 
atoms. One should note that for both sites the Co atom 
is already somewhat below the 1 st Ge layer, having a in- 
coordinate comparable to that of the down-atom for site 
(2) and 0.1 A lower than the down-atom for the site (3). 
As was the case for site (7), these Co atom sites do not 
destroy the reconstruction of the Ge (111)2x1 surface. 

Next, both quasi-stable Co positions were taken as the 
starting positions for additional first principles DFT cal- 
culations with a 4x2 SC and using the same parameters 
as described above. The Co atom was forced to "move" 
into the bulk of the Ge by sequential decrements of its 
z-coordinate (Az — 0.04 A for each geometry relaxation 
step). After each forced sequential decrement of the in- 
coordinate, the position of the Co atom is kept fixed, 
while the Ge atoms are allowed to relax and the total 
energy is determined. By monitoring the total energy of 
the 4x2 SC during subsequent geometry relaxation steps, 
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FIG. 9. (Color online) (a), (b) Calculated quasistable geom- 
etry (4x2 SC) for Co sites (3) and (2) at the Ge(lll)2xl 
surface (see Table [j|. 



we are able to evaluate the potential barrier height that 
needs to be overcome by the Co atom when diffusing ei- 
ther from site (2) or from site (3) to site (7) without de- 
stroying the 2x1 surface reconstruction. First, when the 
Co atom is forced to move only very slightly below the 
Ge surface, the Co atom "bounces" back to its initial po- 
sition. Second, upon a certain minimum translation Az 
along the z-direction (corresponding to a potential bar- 
rier AE), the released Co atom continues to move further 
below the Ge surface to site (7) [see Figs. 10 (a) and (b)]. 



Both Co atom trajectories are visualized schematically 
in Fig. [TO] (a). 

The variation of the total energy E to t of the 4x2 SC 
(again with respect to Etot of the 4x2 SC with the Co 
atom at site (7)) and the variation of the z-coordinate of 
the Co atom during the subsequent geometry relaxation 
steps are presented in Fig. 10 (c) and in Fig. 10 (d) for the 
two considered routes. Three regimes can be discerned. 
Regime I describes the energy gain and z-coordinate vari- 
ations during movement of the Co atom from vacuum to 
either site (2) or site (3). Regime II corresponds to the 
"forced movement" of the Co atom into the bulk of the 
Ge up to the "point of no return" when the potential 
barrier has been overcome. Finally, regime III reflects 
the subsequent relaxation of the Co atom towards its fi- 
nal and stable position at site (7) in the 7-member Ge 
ring. The number of geometry relaxation steps for route 
(2)=K7) and for route (3)=>(7) is 500 and 700, respec- 
tively. For clarity, geometry relaxation intervals where 
the total energy E tot and the z-value remain almost con- 
stant (at the end of the regimes I and III) are cut from 
Figs. 10 (c) and (d). The graphs presented in Fig. 10 (c) 
allow us to determine the potential barrier that the Co 
atom needs to overcome for route (2)=>(7) [(3)=>(7)]: 
AE 2 7 = 0.46 eV [AE 37 = 0.51V], corresponding to a 
change in z-coordinate Azxj = 0.48 A [AZ37 = 0.40 A] in 
Fig. ~ 



10 (d). 



To overcome the surface potential barrier for penetra- 
tion below the Ge surface, a Co atom must have a suf- 
ficiently high (kinetic) energy upon deposition. In our 
experiments Co atoms are evaporated using an e-beam 
evaporator, where the Co mate rial i s heated to a high 
temperature T v around 3000 KPSHSl Atoms leaving an 
e-beam melt generally have a narrow energy distribu- 
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FIG. 10. (Color online) Two different routes are possible 
for the noninvasive embedding of a Co atom from the vac- 
uum side of the Ge surface into the big 7-member ring of the 
Ge(lll)2xl surface reconstruction [also see Fig.[8](a)]: From 
position (2) to (7) and from position (3) to (7). (a) Side and 
(b) top view of the chain-left isomer of the Ge(lll)2xl sur- 
face reconstruction, (c) Relative changes of the total energy 
of the Ge(lll)2xl 4x2 SC (including the Co atom) during 
subsequent geometry relaxation steps for both routes, (d) 
Corresponding variation of the z-coordinate of the Co atom. 



tiorPEI an H the kinetic energy of the evaporated cloud 
of Co atoms, which is induced by the high temperature 
evaporation process can be roughly estimated using the 
equipartition theorem, yielding a mean atomic kinetic 
energy of about 0.38 eV. Recently, however, Asano et 
al. have demonstrated experimentally that the velocity 
of evaporated atoms is typically even higher than the 
maximum velocity suggested by the ideal gas approxi- 
mation.-^ The increased velocity can be accounted for 
by a conversion of electron excitation energy to kinetic 
energy during the adiabatic expansion away from the 
heated material. If the gas flow of evaporated atoms 
cooled sufficiently during the adiabatic expansion, the 
resulting maximum velocity can be estimated as^ 



m 7 - 



1 



RT V , 



(1) 



where 7 is the specific heat ratio C p /C v (C p and C v are 
the specific heat at constant pressure and at constant vol- 
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ume per mole, respectively), R is the gas constant and 
m is the molar mass of the evaporated atom. 7 is 5/3 
for an ideal mono-atomic gas. For T v ~ 3000 K, Eq. (fTl) 
yields a maximum velocity w max = 1460 m/s and hence 
the kinetic energy of the deposited Co atoms may exceed 
even a maximum value of 0.63 eV. This kinetic energy 
allows a Co atom to overcome the surface energy barrier 
AE27 = 0.46 eV or AE 37 = 0.51 eV that is encountered 
when penetrating below the Ge(lll)2xl surface follow- 
ing the route (2)=>(7) or the route (3)=^(7), respectively 
[see Fig. IlO] (c)]. 



E. 



DFT-based modeling of STM topography 
images 



Our theoretical findings, related to the changes in 
surface energy and the potential barrier for penetra- 
tion below the surface, support the idea of a noninva- 
sive embedding of individual Co atoms in subsurface 7- 
member Ge rings. To verify the proposed "embedding 
model" , we investigated the electronic properties of the 
Co/Ge(lll)2xl system by simulating the correspond- 
ing STM topography images using DFT based calcula- 
tions for a wide range of voltages, which allows for a 
direct and detailed comparison between theory and ex- 
periment. For this purpose we investigated the electronic 
structure of the Ge(lll)2xl surface for each of the pos- 
sible Co atom locations described in Section |IVB[ for 
both chain-left and chain-right isomers. The calculated 
quasi-stable geometries of the relaxed 8x4 SCs (9 and 8 
possible geometries for the chain-left and chain-right iso- 
mer, respectively) were transferred to a larger 14x4 SC 
(size is 55.82 A x 27.63 A) for which we calculated the 
electronic properties in detail. 

In order to construct STM constant current topogra- 
phy images based on the calculated electronic structure 
of the Co/Ge(lll)2xl system, the decay of the electron 
wave functions from the surface into the vacuum needs to 
be taken into account. Within the Tersoff-Hamann the- 
orySMSl the dependence of the tunneling current / on the 
applied tunneling voltage Vt between an STM tip and a 
surface is given by 



27re y^ 



/(£„)[! - f(E v + eV t ] I M^ | 2 $(£„ - E v ), 



(2) 

where f(E) is the Fermi function, M^ v is the tunneling 
matrix element between electronic states ip^ of the tip 
and electronic states ip v of the surface, and E^ (E v ) is the 
energy of the state ip^ (ip v ) in the absence of tunneling. 
When we assume localized wave functions ipa for the tip, 
M^ will vary proportional to the amplitude of ip v at 
the position r 0, which corresponds to the center of the 
sphere that is used to approximate the tip apex. At low 
temperatures and for small tunneling voltages V t Eq. p]) 
reduces to 



I<xJ2\M?o) 



8{E^ - E F ). 



(3) 



From Eq. p| it follows that the tunneling current / is 
proportional to the surface LDOS that is probed at posi- 
tion r of the tip, integrated over an energy range from 
Ep to E-p + eVt- For constant tunneling current I = It 
the STM tip essentially follows a contour of constant sur- 
face LDOS. However, because the surface wavefunctions 
decay exponentially into the vacuum region, numerical 
evaluation of ijJ v (ro) for tip-surface distances of the or- 
der of several angstroms becomes a significant problem 
for DFT calculations!^ For this reason STM simulations 
are often restricted to (the vicinity of) the surface, which 
may yield incorrect results. To tackle this problem, we 
have used the 2D Fourier transform of the wavefunctions 
ijj v (r) in combination with spatial extrapolation tech- 
niques^ to evaluate the surface wave function ip v (x, y, z) 
in the vacuum region up to z = 7 A above the surface. 

This way we calculated the STM topography images 
for all available Co atom sites and for both the chain-left 
and chain-right isomers (for a 14x4 SC) within an energy 
range between — 1.5eV and +1.5 eV and at distances up 
to 7 A above the surface. Perfect agreement between the- 
ory and experiment for the whole energy range can be 
achieved only for a Co atom located at site (7) for the 
chain-left isomer geometry. For site (7) we then calcu- 
lated the electronic properties also for a 26x4 SC (size 
is 103.67 A x 27.63 A, consisting of 1873 atoms) and a 
23x5 SC (size is 91.74A x 34.54A, consisting of 2071 
atoms). CG geometry optimization for the chain-left(2) 
isomer with the Co atom located at site (7) within a 
9x3 (35.89 A x 20.72 A) SC was pe rform ed, similar to 
the calculations described in Section IIVBI We used the 



chain-left (1) and chain-left (2) isomers (see Section IV A) 



for the 26x4 SC and 23x5 SC, respectively, to keep the 
Co atom in the center of the SC. Calculations for the 
26x4 chain-left (1) SC and the 23x5 chain-left (2) SC 
yield identical results. 

For our simulations of the STM topography images we 
have to rely on experimental z(Vt) spectra measured on 
the Ge(lll)2x 1 surface in order to take into account the 
dependence of the height z on the tunneling voltage V t 
in our calculations. The experimental z(V t ) dependence 
with an initial height addition of 3 M^was used to deter- 
mine the height above the Ge(lll)2xl surface at which 
simulated STM images are calculated. For low voltages, 
i.e. for energies close to Ep , z(V t ) ~ 3.7A, while for high 
voltages above IV, z(V t ) ~ 6 A. 

In Fig. 11 we present a series of experimental (inner 
columns) and calculated (outer columns) STM topogra- 
phy images for the filled (two right columns) and empty 
(two left columns) states regime within a wide range of 
tunneling voltages. Calculated STM topography images 
are obtained for a 23x5 SC with chain-left (2) isomer ge- 
ometry with the Co atom at site (7). Experimental STM 
topography images are all recorded at the same location. 
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FIG. 11. (Color online) 9.2 x 3.5 nm 2 experimental empty and filled states STM topography images (inner columns), together 
with the corresponding calculated DFT-based STM topography images (outer columns) of a single Co atom located in a 
subsurface 7-member Ge ring, in between the 3 rd and 4 th atomic layer underneath the Ge(lll)2xl surface. The tunneling 
voltage Vt is indicated for each image. 



The Co atom is well separated from other Co atoms, im- 
lying there is no influence from neighboring Co atoms 
(see Section III). 



As can be seen in Fig. [TTJ correspondence between the 
Co related features in the calculated and the experimen- 
tal STM topography images is striking for the whole in- 
vestigated voltage range. Concerning the precise tun- 
neling voltage V t at which optimum correspondence is 
observed between theory and experiment there is a mi- 
nor mismatch. This mismatch exhibits a non-linear de- 
pendence on the applied tunneling voltage for both the 
filled and empty states regime. At low tunneling volt- 
ages the difference in voltage is around 0.1 V for both 
the empty and filled states regime (see the first row of 
images in Fig. 111. For higher tunneling voltages, the 



difference increases to around 0.37 V and 0.23 V for the 
filled and empty states regime, respectively. Upon more 



careful comparison, it can be seen that the difference 
in voltage mainly affects the surrounding Ge(lll)2xl 
surface and not the Co atom itself. Indeed, maximum 
contrast related to the Co protrusion appears around 
V t = 0.9 ± 0.1 V in both the experimental [Fig. 11 and 
Fig. [6] (a4)] and the calculated STM topography images. 
Apart from the rather small difference in voltage, there 
is a very good agreement between the theoretically and 
experimentally observed electronic features for both the 
filled and empty states regime. 

The results presented in Fig. [TTJ confirm that the cal- 
culated images of the Co/Ge(lll) system exhibit all typ- 
ical features that were obser ved in the voltage dependent 
STM investigation (Section |!!lC| ): (i) The Co induced 
strongly perturbed area comprises ±2 SUCs on the up- 
per 7r-bonded chain row at moderate and high energies 
in the empty states regime; (ii) the perturbation has a 
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FIG. 12. (Color online) Calculated empty state STM topogra- 
phy image of a Co atom located at site (7) in the Ge(lll)2xl 
surface. The empty state energy (0.57 V)is close to the energy 
for optimum calculated contrast of the zigzag structure of the 
upper 7r-bonded chains. 




mirror symmetry axis along the [211] direction; (iii) the 
perturbation exhibits a clear asymmetry with respect to 
the [Oil] direction; (iv) the Ge(lll)2xl surface exhibits 
a zigzag structure only at a specific tunneling voltage and 
this zigzag structure is perturbed near the Co atom. 

In the filled states regime the calculated ID Co induced 
perturbation along the 7r-bonded chain row exhibits the 
highest contrast in the —0.25 eVto — 0.35 eV voltage 
range, whereas experimentally the highest contrast oc- 
curs around eV t — —0.65 ± 0.05 eV [see Fig. [6] (b3) and 
Fig." 
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At these energies the Co induced perturbed 
atomic corrugation along the upper 7r-bonded chain row 
in the calculated and experimental STM topography im- 
age match very well. At energies below —1.0 eV the Co 
atom appears as an atomic size vacancy in the upper 
7r-bonded chain row for both the calculated and experi- 
mental STM topography. 

In the calculated images the zigzag structure of the 
Co containing upper 7r-bonded chain rows is observed at 
empty state energy around 0.57eV (see Fig. 12 1. On 



the other hand, this zigzag structure appears around 
Kigzag = 0-8 5 eV in the experimental STM images (see 
Section III B I . The determination of the angle a in Fig. 12 
allows for an easy comparison to Fig. H (d). It can be 
seen that, apart from a small energy mismatch, the Co 
induced perturbation of the calculated zigzag structure 
again perfectly matches the experimental observation. 

The energy mismatch between the calculated and ex- 
perimental STM results can be related to the doping of 
the investigated Ge crystals (p-type Ge crystals with a 
low dopant concentration are used in this work), which 
is not included in the DFT modeling. Since the surface 
and bulk bands shift in energy depending on the type of 
doping and on the doping level, it can be expected that 
the Co induced perturbations shift in energy as well. On 
the other hand, whereas deposition of 0.032 ML of Co did 
not lead to changes of the Ge(lll)2xl electronic struc- 
tures in the experiments, the incorporation of the Co 
atom in the DFT model may induce a "doping" effect of 
the Co/Ge(lll)2xl system due to the finite size of the 
SC. For the 23x5 SC, the Co/Ge ratio is 1/2071, which 
corresponds to a heavily doped Ge crystal. As already 
mentioned in Section |HIB| the precise value of the tun- 



FIG. 13. (Color online) (a), (b) 14.0 x 21 nm 2 experimen- 
tal filled states STM images of 6 individual Co atoms for 
n-type Ge(lll)2xl (resistance pbuik — llficm). The tun- 
neling voltage Vt is indicated for each image. It is fixed at 
100 pA for (a) and at 500 pA for (b). Blue cross markers in- 
dicate identical locations in (a) and (b). 



neling voltage 14ig Za g is found to depend on the semicon- 
ductor type and on the doping concentration. E.g., heav- 
ily doped n-Ge(lll)2xl surfaces (phosphor doping level 
rip = 1 x 10 19 cm~ 3 , surface preparation as described in 
Section [n]) are found to have V Z i gzag — 0.55 V, which is 
in good agreement with the calculated tunneling voltage 
of 0.57V. 

Finally, we investigated the influence of the doping 
level on the appearance of the clean Ge(lll)2x 1 surface 
and the Co/Ge(lll)2x 1 system in experimental STM to- 
pography images. As indicated above, it can be expected 
that the Co induced perturbations shift in energy, since 
the surface and bulk bands shift in energy depending on 
the type of doping and the doping level. In Fig. [13] we 
present two experimental filled states STM topography 
images of a low doped ra-type Ge(lll)2xl surface. The 
Ge crystal was doped by P at a doping level np = 0.5 
to 1.0 x 10 15 cm~ 3 (pbuik — UCl cm). Surface prepa- 
ration and Co deposition are performed as described 
in Section [TTl At low temperatures (T samp \ c ~ 4.5 K), 
these samples exhibit a detectable tunneling current in 
the filled states regime only at tunneling voltages below 
—0.8 V. Interestingly, the above described Co induced 
perturbations appear at significantly different energies 
for this low doped sample. The features observed in 
Fig. 13 (a) for the n-type Ge(lll)2xl surface are ob- 
served around V t — —0.7 V for the p-type Ge(lll)2xl 



surface (see Fig. 11 ), corresponding to an energy shift of 
about 0.7 to O.i 

The energy 



ieV. A similar energy shift can be inferred 
by comparison of Fig. 13 (b) and Fig. 
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shifts between the STM experiments and the DFT cal- 
culations in Fig. [TTJmay therefore be attributed to doping 
effects. Alternatively, the observed energy shifts may also 
be related to the intrinsic deficiency of LDA with respect 
to the quantitative determination of band gap values of 
semiconductor materials. 
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FIG. 14. (Color online) Chemical bonding in terms of the 
covalent bond energy (Ecov) f° r Ge-Ge bulk and Co-Ge in- 
teractions. 



F. Co/Ge(lll)2xl electronic properties and 
Co— Ge bonding characteristics 

In this section we discuss the chemical bonds that are 
formed between the Co atom at site (7) and its surround- 
ing Ge atoms [labeled in Fig. 15 (bl)] underneath the 



Ge(lll)2xl surface. The amount of possible chemical 
bonds can be roughly estimated by relying on a simple 
analysis of the Co-Ge bond lengths. For a single bond 
between a Co and a Ge atom the bond length can be 
estimated as the sum of the Co and Ge covalent radii, 
which is r c = 2.38A. The calculated distances between 
the Co atom and the neighboring Ge atoms for the 8x4 
SC equilibrium geometry (see subsection IV B ) are listed 
in Table plj Among the listed Ge atoms, Ge atoms (2) 
to (5) are most likely to form a covalent bond with the 
Co atom. The distances between the Co atom and these 
Ge atoms are, however, slightly larger (1-2%) than r c , 
which can be attributed to the employed LDA (an over- 
estimation of the bond lengths by a few percent is typical 
for LDA). On the other hand, the calculated Co-Ge(l) 
distance is significantly larger (around 13%) than r c , im- 
plying that the formation of a bond between the Co atom 
and the Ge(l) wp-atom can be excluded. 

An estimate of the atomic charges, electron transfer 
and covalent interactions between the Co atom at site (7) 
and the neighboring Ge atoms [labeled in Fig. 15 (bl), 
also see Table In] can be obtained more qualitatively 



TABLE II. Calculated distances between the Co atom at site 
(7) and the neighboring Ge atoms [labeled in Fig.[l5](bl)] for 
the 8x4 SC equilibrium geometry. 



Co neighboring atoms 


Co-Ge distance 


Co-Ge(l) 


2.687 A 


Co-Ge(2) 


2.404 A 


Co-Ge(3) 


2.404 A 


Co-Ge(4) 


2.415 A 


Co-Ge(5) 


2.243 A 



from the calculated number of electrons inside the atomic 
spheres and from the Mulliken populations.^ The Mul- 
liken electron orbital overlap populations (calculated for 
the 23x5 SC) of Ge atoms (1) to (5) both with (per- 
turbed 7r-bonded chain row of the Co/Ge(lll)2xl sys- 
tem) and without (unperturbed 7r-bonded chain row of 
the Ge(lll)2xl surface) the Co atom are given in Ta- 
bic [III] (A) and in Table [LlITB), respectively. The (un- 
perturbed) bulk atoms Ge(4 ) and Ge(5') have four Ge- 
Ge covalent bonds, with an electron orbital overlap pop- 
ulation around 0.40 ± 0.01 e/bond (sp 3 hybridization). 
On the other hand, up-atom Ge(l') and down-atoms 
Ge(2',3') have only three covalent bonds (sp 2 hybridiza- 
tion): A fourth bond is absent or at least strongly re- 
duced (Mulliken overlap populations are less than 22% 
of the Ge-Ge bulk covalent bond). Co-induced changes 
in the Mulliken overlap population of the neighboring Ge 
atoms and the atom population can be traced by com- 
parison of Table[Hl](A) to Table|Hi](B). The calculations 
reveal that no significant charge redistribution between 
the Co atom and the surface and bulk Ge atoms occurs, 
i.e. the electron population of the Co atom remains close 
to that of a neutral Co atom (9.043 e). The Co atom 
exhibits four covalent bonds with Mulliken overlap pop- 
ulations that are 56% to 60% of the Ge-Ge bulk covalent 
bond. There exists a weak interaction with the Ge(l) up- 
atom as well (see Table III (A)). Because of the very low 
overlap population (30% of Ge-Ge bulk bond) and the 
very weak charge transfer between Co and Ge(l), how- 
ever, the interaction between the Co and Ge(l) up-akom. 
should not be considered as a fifth covalent bond. 

The energy of the covalent Co-Ge bonds can be esti- 
mated as the total energy difference between a Co atom 
that is "bonded" [Co atom located at site (7)] or "not 
bonded" to the Ge lattice. Calculations for the "not 



TABLE III. (A) Mulliken overlap population and electron 
population of the perturbed Ge atoms neighboring the em- 
bedded Co atom at site (7) and (B) of the unperturbed Ge 
atoms in the absence of the Co atom (ideal Ge( lll) 2xl sur- 
face). Ge atoms are numbered according to Fig. 15 (bl) and 
Table [III 



A. Co/Ge(lll)2xl 
Bond (atom-atom) 

Ge-Ge 

Ge-Ge 

Ge-Ge 

Ge-Ge 

Co-Ge 
Electron pop. 



Ge(l) Ge(2) Ge(3) Ge(4) Ge(5) 



0.368 
0.394 
0.396 

0.120 
4.205 



0.444 
0.394 
0.409 
0.067 
0.241 
3.891 



0.396 
0.445 
0.408 
0.065 
0.239 
3.889 



0.368 
0.401 
0.228 
0.401 
0.224 
3.943 



0.228 
0.356 
0.401 
0.353 
0.231 
3.884 



B. Ideal (2x1) 
Bond (atom-atom) 



Ge(l') Ge(2') Ge(3') Ge(4') Ge(5') 



Ge-Ge 
Ge-Ge 
Ge-Ge 
Ge-Ge 
Electron pop. 



0.412 
0.476 

0.475 

4.172 



0.475 
0.476 
0.435 
0.089 

3.882 



0.475 
0.456 
0.435 
0.089 
3.859 



0.412 
0.402 
0.395 
0.402 
3.982 



0.395 
0.406 
0.405 
0.406 
3.961 
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0.077 e/A 3 



0.059 e/A 



0.054 e/A 3 



0.050 e/A 3 



FIG. 15. (Color online) 3D isosurface charge density p(x, y, z) plots of the Co/Ge(lll)2xl system (chain-left isomer) with the 
Co atom located inside the big 7-member Ge ring below the surface. Viewpoint is along the [0111 direction for (al) to (a4) 
and along the direction indicated by the arrow A in (al) for (bl) to (b4). The charge density p (e/A 3 ) isovalues for (al), (a2), 
(a3) and (a4) are the same as the ones for (bl), (b2), (b3) and (b4), respectively, and are indicated at the bottom of the latter 
parts of the figure. 



bonded" case were performed for a 4x2 SC with the Co 
atom placed 4A above the Ge surface. The total energy 
difference is found to be 8.4 ± 0.3 eV. Considering four 
Co-Ge bonds as discussed above, the (average) energy of 
a single Co-Ge bond is hence -Eco-Gc = 2.1 eV, which is 
significantly lower than the energy of the covalent Ge-Ge 
bond (£ Gc _ Gc = 3.71 eV). 

The analysis of the bond lengths and the Mullikcn 
overlap populations are indicative of rather weak Co-Ge 
bonds, which can also be concluded from our experimen- 
tal observation that Co atoms diffuse even at lower tem- 
peratures!^] More precisely, in spite of the low sample 
temperature (T samp i c < 80 K) during Co deposition, the 
majority (around 87%) of the Co atoms diffuses along ir- 
bonded chain rows to surface/subsurface defects (includ- 
ing DBs, MASs, and subsurface Ga impurities), which 
can be related to the weak Co-Ge bonds. The remain- 
ing 13% can be retrieved as individual (subsurface) Co 
atoms at the cold Ge surface and they diffuse as well to 
surface/subsurface defects after warming the substrate 
up to room temperature. For a more detailed discussion 
we refer the reader to Ref. p^. 

The chemical bonding mechanism between the Co and 
Ge atoms involved in the DFT electronic structure calcu- 
lations can be investigated in more detail by evaluating 
the crystal orbital overlap population/Hamiltonian pop- 
ulation (COOP/COHP)P125| We have used an a i terna . 

tive COOP/COHP based approach that allows to calcu- 
late the relevant physical quantities independent of the 
choice of zero of the potential b y re lying on the so-called 
"covalent bond energy" (-Ecov)- COOP and -Ecov cal- 



culations are known to yield similar results, while the 
COOP method generally overestimates the magnitude of 
the anti-bonding states when defined within a plane- wave 
basis set .^S Figure 14 presents our Eqov calculation for 
the Co-Ge and Ge-Ge bulk interactions in an energy 
range of 12 eV around the Fermi level. Note that the 
Ecow values (y-ordinate) are plotted without any units 
and can only be interpreted qualitatively. Negative, pos- 
itive, and zero values of -Ecov correspond to bonding, 
anti-bonding, and non-bonding interactions, respectively. 
The Ecov spectra confirm the stability of the Co/Ge sys- 
tem. Strong bonding interactions exist for Ge-Ge from 
the bottom of the VB up to the Fermi level, while for the 
Co-Ge a strong bonding interaction is found only up to 
around —0.3 eV. Between —0.3 eV and the Fermi level, 
Co-Ge anti-bonding interactions occur, which may ex- 
plain the experimentally observed thermal instability of 
the Co/Ge system (the covalent Co-Ge bonds are weaker 
than the Ge-Ge bulk bonds). 

Finally, we calculated charge electron density maps to 
visualize the Co-Ge bonds. In Fig.[l5]we present isosur- 
face maps of the spatial electron charge density p(x, y, z) 
for isosurface values ranging from 0.077 e/A 3 down to 
0.050 e/A 3 . High electron densities related to the Ge- 
Ge covalent bonds gradually appear above 0.059 e/A 3 
in Figs. fl5] (a2,b2) to (a4,b4). Moreover, zones of high 
electron localization between the Co atom and the sur- 
rounding Ge atoms can be observed. A zone of high 
electron localization exists between the Co atom and the 
Ge(5) atom in Fig. [15] (a2) (indicated by red arrows). 
Two additional symmetrical zones of high electron local- 



ization can be observed between the Co atom and the 
Ge(2) and Ge(3) down-atoms in Fig. Kb\ (b2) (indicated 
by the blue arrows). A fourth Co-Ge bond can be re- 
lated to the high electron localization zone between the 
Co atom and the Ge(4) atom in Fig. 15 (a3) (indicated by 
the green arrows). As could be expected from our above 
analysis, there occurs no high electron localization zone 
between the Co atom and the up-atom Ge(l), which ad- 
ditionally confirms that both atoms are not bonded. The 
findings related to the bond lengths, the Mulliken over- 
lap populations and the -Ecov calculations are hence in 
agreement with the isosurface charge density maps. We 
therefore conclude that the Co atom forms a bond with 
four neighboring Ge atoms, corresponding to a Co 4+ va- 
lence state and a 3d 5 electron configuration. It is known 
that Co has a wide range of valence states due to its var- 
ious spin configurations, implying Co 4+ . i.e. 3d 5 ions can 



exist in several spin configurations ' 98 ^ 99 ^ An extension of 
our DFT model for the embedding of a Co atom in the 
Ge lattice by including spin-dependent calculations will 
be a topic of future research. 



CONCLUSION 
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surface and clear ID confinement along 7r-bonded chains. 
Relying on DFT based calculations we demonstrated that 
the energetically most favorable position of a Co atom 
is attained by penetration into the Ge(lll)2xl surface. 
The Co atom occupies a quasi-stationary position within 
the big 7-member Ge ring of the Ge(lll)2xl reconstruc- 
tion in between the 3 rd and 4 th atomic layer beneath the 
surface. The embedded Co atom induces an electronic 
asymmetry of the 7r-bonded chain with respect to the 
[Oil] direction, which allows us to determine that the 
Ge(ll 1)2x1 surface reconstruction of all investigated Ge 
samples consists of domains with chain-left geometry ex- 
clusively. Calculated STM images based on our DFT ap- 
proach match very well the experimental STM images of 
the Co/Ge(lll)2x 1 system within the investigated range 
of tunneling voltages. Finally, DFT based calculations of 
the Co-Ge bond strength reveal the formation of four 
covalent bonds, corresponding to a Co 4+ valence state 
and a 3d electron configuration. Our findings open in- 
teresting perspectives for investigations of subsurface ID 
(nanowires) and 2D (islands) nanostructures that are ex- 
pected to form at a higher coverage of Co on Ge(lll). 
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